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Abstract: The quark-gluon vertex in Landau gauge is studied in the quenched approx- 
imation using the Sheikholeslami-Wohlert (SW) fermion action with mean-field improve- 
ment coefficients in the action and for the quark fields. We see that the form factor that 
includes the running coupling is substantially enhanced in the infrared, over and above 
the enhancement arising from the infrared suppression of the quark propagator alone. We 
define two different momentum subtraction renormalisation schemes — MOM (asymmet- 
ric) and MOM (symmetric) — and determine the running coupling in both schemes. We 
find A— L = 300jli8o ± 55 ± 30 MeV from the asymmetric scheme. This is somewhat 
higher than other determinations of this quantity, but the uncertainties — both statistical 
and systematic — are large. In the symmetric scheme, statistical noise prevents us from 
obtaining a meaningful estimate for A^jg. 
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1. Introduction 

The quark-gluon vertex plays an important role in many applications of QCD. QCD vertex 
functions may be used to define momentum subtraction (MOM) renormalisation schemes 
These, it is argued, are more 'physical' than the minimal subtraction (MS or MS) 
schemes, since the latter can only be defined in a perturbative context, while the former are 
independent of the regularisation method and give a better guidance to the appropriate 
renormalisation scale for a particular problem. Recently, a complete determination of 
the quark-gluon vertex to one-loop order was performed In an asymmetric momentum 
subtraction scheme, it has recently been computed to three-loop order , while a numerical 
computation to two-loop order has been performed in a symmetric scheme ||. 
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The quark-gluon vertex from the lattice may thus yield a direct determination of the 
QCD running coupling a s , to complement other methods for determining this quantity Q. 
For a review of experimental and theoretical determinations of a s , see J?]]. 

In addition to reproducing the ultraviolet, perturbative behaviour, and thus deter- 
mining the intrinsic QCD scale Aj^g, the quark-gluon vertex can also be used to probe 
the infrared behaviour of the running coupling. The hypothesis that the QCD coupling 
approaches a constant in the infrared has long been popular on phenomenological grounds, 
and has also received some theoretical support from 'optimised' perturbation theory (8|. 
On the other hand, a more recent proposal for reorganising perturbation theory H gives 
a running coupling that vanishes in the infrared. Lattice QCD may in principle assist in 
resolving this issue. 

The nonperturbative quark-gluon vertex also enters into the Dyson-Schwinger equa- 
tions (DSEs) 11, 12 1 , which are the QCD field equations. In particular, in Minkowski 



space the DSE for the renormalised quark propagator S(p) is 



Z 2 (j* - Z rn m) + i^Z 1F g 2 J 



d 4 k 



7tl S(k)D^(k-p)T u (p,k-p) 



XI) 



where m and g are the renormalised mass and coupling respectively. The unknown quan- 
tities here are the nonperturbative gluon propagator D^ u (q) and the quark-gluon vertex 
(p, q) [see figure [l]] . Z2 and Z m are the quark field and mass renormalisation constants 
respectively. It will be convenient to also introduce A u (p,q) = —igT u (p,q). 

Since these are gauge dependent quantities, 
it is to be expected that the confinement pic- 
ture, and the relative importance of the various 
factors, will vary between different gauges. We 
will here be working in the (minimal) Landau 
gauge, where over the past few years substan- 
tial progress has been made in our understanding 
of the gluon propagator from lattice simulations 

m 



fi] a 

a 



14, 15, uM as well as analytical studies fll 



k = p + q 



Figure 1: The quark-gluon vertex. 



The quark self-energy directly exhibits con- 
finement and dynamical chiral symmetry break- 
ing, and is an important input for phenomeno- 
logical models of hadron physics [[U]] . It has also 

recently been evaluated in Landau gauge on the lattice [[0| [0|, 20,21]. 

However, the quark-gluon vertex remains largely unknown, and the validity of the 
usual ansatze untested. In Landau gauge, there are indications that it must contain non- 
trivial structure in the infrared. The lattice (infrared suppressed) gluon propagator together 
with a bare or QED-like vertex fails to yield solutions to ( |l.l|) that exhibit an appropriate 



degree of dynamical chiral symmetry breaking [22|. But there are also strong indications 
that the ghost propagator in Landau gauge is strongly enhanced, both from lattice simula- 
tions 24 1 and analytical studies [17|. The ghost self-energy enters into the quark-gluon 
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vertex through the Slavnov-Taylor identity, 

gT M (p, q) = G(q 2 ) [(1 - B(q,p + q^S^ip) - S^ip + q)(l - B(q,p + q))] , (1.2) 

where G(q 2 ) is the ghost renormalisation function and B a (q,k) = t a B(q,k) is the ghost- 
quark scattering kernel, which is given by the diagram in figure ||. It appears that modelling 
this into the quark DSE does give solutions exhibiting chiral symmetry breaking and quark 
confinement [25|] . 

A nonperturbative determination of the quark-gluon vertex will therefore give us fur- 
ther insight into the mechanisms of confinement and chiral symmetry breaking, as well as 
casting light on the transition between the perturbative and nonperturbative regimes of 
QCD. 

The starting point for analytical studies 
of the quark-gluon vertex is the QED vertex, 
which gives the 'abelian' contribution to the 
QCD vertex. In the abelian case, the Slavnov- 
Taylor identity implies that the vertex can be 
written entirely function of the nonper- 

turbative quark propagator up to a transverse 
term, as shown by Ball and Chiu (2(| . A kine- 
matical basis, along with a one-loop determi- 





q; a 



Figure 2: The ghost-quark scattering kernel 
B a . 



nation of all the components, is given in [26 

H- 

The three-gluon vertex has been the sub- 
ject of detailed study in recent years |28|, |2t| 
An important result of these studies is the discovery of substantial power corrections to 
the running coupling, which remain up to scales of 7-10 GeV, and originate from the {AA) 
condensate appearing in the Landau gauge OPE. It has been conjectured that this conden- 
sate is largely due to instanton effects |30j j. In the quark-gluon vertex the situation is in 
some senses more straightforward. Power corrections due to the quark mass appear already 
in the one-loop perturbative running coupling, and these are expected to be substantially 
enhanced by the chiral condensate. These are phenomena that will appear in any gauge. 
Any correction due to a covariant-gauge {AA) condensate will come in addition to this. 

This paper builds on earlier results that were presented in |[l] . We will be focusing on 
the form factor multiplying 7„, which contains the running coupling. We investigate the 
infrared behaviour of this form factor, to study the hypothesis of infrared enhancement, as 
well as the ultraviolet behaviour, attempting to determine the perturbative running. 

In section |2| we define our notation and the quantities involved. In section |3] we define 
two momentum subtraction schemes based on the quark-gluon vertex, which we call MOM 
and MOM, and explain how the running coupling may be extracted in these schemes. The 
parameters of our simulations are given in section f|. In section [5] we determine the quark 
and gluon field renormalisation constants Z2 and Z3. Our results for the Ai form factor and 
the MOM running coupling are given in section 6T, while the results for the MOM scheme 
are given in section |6.2j . In section [j] we translate the MOM scheme results to the MS 
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scheme. Finally, in section |8| we summarise our conclusions and discuss the prospects for 
further work. The appendices contain a discussion of the full tensor decomposition of the 
vertex; the lattice tree-level expressions and our method for removing the dominant (tree- 
level) lattice artefacts; and the full one-loop continuum expressions for the form factors we 
consider. 

2. Definitions and principles 

With the exception of section [7|, where the perturbative expressions are given in Minkowski 
space, we will be working throughout in euclidean space, with a positive metric, such that 
A 2 > for any spacelike vector A. The commutation relations for the Dirac matrices are 
the usual ones, 

{ 7m , 7 ,} = 2V; 7J = 7m- (2-1) 
The (7^-matrices are defined by 

<V = 2 (2-2) 

The generators t a of the Lie algebra have the conventional normalisation, Tr(t a t b ) = ^S ab . 

We can define the configuration space quark-gluon vertex function (see figure |T|) on 
the lattice as 

V2(x, y ,z)% = (<(x)^(z)^(y)) = ((s%(x,z)A^y))) . (2.3) 

Here, (• • •) denotes averaging over all fermion and gauge field configurations, while ((• • •)) 
denotes averaging over gauge field configurations only. Fourier transforming this and invok- 
ing translational invariance gives us the full (unamputated) momentum space bare vertex 
function V£(p,q): 

J2e-^-y-^(^x)A%y)^z^ 

x,y,z 

= e -*(P-*+*)- e-^™) (<(*)^(y)^(0)) 

z x,y 

= V5{p-k + q) ((s%(p;U)A^q))) = V5{p-k + q) V^(p,q)%, (2.4) 

where V is the lattice volume. The proper (one-particle irreducible) bare vertex can be 
obtained by amputating the external quark and gluon legs from the full vertex V^: 

kf*\p,q) = S(p)- 1 V-(p, q )S(p + q)- 1 D( q )^ . (2.5) 

The only possible dependence this can have on the group coordinates a, i,j is proportional 
to the generator £?.. We can therefore consider only A^{p,q), defined by 

= *tf( A A«W» = -*$btyMa0 • (2.6) 
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S(p) = ((S(p;U)}} = ((Ylx e ipx S(x,0;U))) is the momentum-space quark propagator, 
while D{q) is the gluon propagator, which in the infinite-volume limit takes the form 

D*(q) = 5 ab D,M) = 5 ab - ^) D(q 2 ) + 6« b ^± . (2.7) 

In the Landau gauge (£=0), D{q 2 ) can be determined for q ^ by 

D(q 2 ) = — -1 -Vrj?). (2. 

W ' 3(N 2 - 1) ^ wW 



As long as q is not too close to zero, this form remains valid also for a finite volume. In 
general, a finite volume will induce an effective 'mass' m ~ 1/L, which on an asymmetric 
lattice may also be direction-dependent — so the tensor structure ( |2.7| ) must be replaced 

by @ 

DfM = 5 ah (V " ^) D(q 2 ) + ^e^y - S a % u (q)D(q) + , (2.9) 



where the functions f,g,h,h' are such that for sufficiently large q, ( |2.9[ ) approaches the 
infinite-volume form, but both / and g remain non-zero as q — > 0. The Landau-gauge 
expression (2.8) must for the smallest momentum values — and in particular for q = on 
any volume — be replaced by 

D[q2) = W^T) £ DZiQ)/ £ 55 r (g )(4-D £ • (2 - 10) 

In the infinite- volume limit, 7^(0) — > 5a V since the Landau gauge condition places no 
restriction on the zero- modes P4] , so T(0) — ► 4. In [13] it was found that an asymmetric 
finite volume may induce large distortions to this form, and in general the components 
must be determined numerically. However, T^(0) will always remain diagonal. 

Since the gluon propagator in Landau gauge for q 7^ becomes proportional to the 
transverse projector P^ u (q) = — Q^Iu/q 2 or its lattice equivalent, D~ l is undefined and 
we cannot use eq. (|2.5|). Instead, we rewrite (2.5) as follows, 



D^(q)Al(p,q) = P^(q)D(q 2 )A a u (p,q) = S(p)- 1 V«(p,q)S(p + q)- 1 , (2.11) 

from which we can obtain the transverse-projected vertex, 

A a /(p,q) =P„u(q)K M (P,<l) =S{p)- 1 V*{p 1 q)S(p + q)- l D(q 2 )- 1 . (2.12) 

The quantities calculated on the lattice are always functions of the bare (unrenor- 
malised) fields tp , A° and the bare coupling g$ . The relation between renormalised and 
bare quantities is given by 

< = Z 2 1/2 V>; / = Z 2 1/ V; Al = zl /2 A^- g = Z g g; £ = Z 3 £, (2.13) 
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where Z2,Z 3 ,Z g are the quark, gluon and vertex (coupling) renormalisation constants 
respectively, and are functions of the regularisation parameter a and the renormalisation 
scale [i. From ( |2.13| ) it follows that 

S b ^(p;a) = Z 2 ^;a)S(p^); (2.14) 
D h ^(q 2 ;a)=Z 3 (^a)D(q 2 ;^. (2.15) 

The renormalised vertex is related to the bare vertex according to 

A£ are (p, q; a) = a)A>, q; /x) . (2.16) 

Gauge invariance requires that Z\p = Z g Z2Z\j 2 , and so we may write 

A h ™ e (p,q;a) = Z^MZ^MZ-^MA^q-,^), (2.17) 

meaning that only the quark and gluon fields, along with the running coupling, are indepen- 
dently renormalised. For the sake of brevity, we will from here on no longer explitly label 
bare quantities as such. Z2 and Z3 may be determined by imposing momentum subtrac- 
tion (MOM) renormalisation conditions on the quark and gluon propagator respectively, 
demanding that they take on their tree-level forms at the renormalisation scale 

S(p;v) | 2 2 = . * -, s ; (2.18) 

D{q 2 ^) \ q2=fl2 = j 2 . (2.19) 

The renormalisation of the coupling will be discussed in section pL 

The Lorentz structure of the vertex in the continuum consists of 12 independent vectors 
and can be written in terms of a 'Slavnov-Taylor' (non-transverse) and purely transverse 
part in terms of vectors hi,T% and scalar functions Aj,Tv 

A^p,q)=Al ST \p,q)+A^(p,q) 

4 8 (2.20) 

= -i9^2 ^i(p 2 ,q 2 ,k 2 )L iifl (p,q) - ig^2Ti(p 2 ,q 2 ,k 2 )T iifl (p,q) . 

i=l i=l 

The full expressions for all the vectors Lj and Tj are given in appendix [A|. In this paper, we 
will only study the part of the vertex proportional to 7^, which in the specific kinematics 
we will be employing is given by the three vectors 

Li ttl = 7 M ; L 2 ^ = -(2 /$+ 4)(2p + g)„ ; T 3 ^ =^ - q 2 7fl . (2.21) 

Because of the Slavnov-Taylor identity ( |1.2|) , the scalar functions Xi(p 2 ,q 2 ,k 2 ) in ( 2.2C| ) 
may be expressed in terms of the quark propagator, ghost propagator and ghost-quark 
scattering kernel. In QED for instance, as a result of the Ward-Takahashi identity, Ai is 
given uniquely in terms of the fermion propag ator S^ip) = i fA{p 2 ) + B(p 2 ) 1 by p§ 

A? ED (pW, k 2 ) = \ (A(p 2 ) + A(k 2 )) . (2.22) 



Note that the opposite conventions for the B-function are often used in Minkowski space: our B 
corresponds to — (3 in H Q 
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In QCD, for q = 0, the equivalent of this is [33] 



Xx(p 2 ,0, P 2 ) = G(0)[a( P 2 )xo(p 2 ,0, P 2 ) + B(p 2 )(xi(p 2 ,0,p 2 ) + X2(p 2 ,0,p 2 )) 

-2p 2 A{p 2 )xz{p 2 ^p 2 )] , (2.23) 

where Xi are the fo rm factors of the ghost-quark scattering kernel given in ||. At tree 
level, xo = l,Xi,2,3 = 0. 

As already mentioned, in Landau gauge, for q 2 ^ 0, we can only determine the trans- 
verse part of the vertex from the lattice. The transverse projection of L\ is 



P^{q)Li, u (p,q) = ~^T 3:tM (p,q) . 

This gives rise to the modified form factor A^, 

A'i = Ai - q 2 r 3 , 

which will be useful when studying the transverse-projected vertex. 



(2.24) 



(2.25) 



3. Definition of the MOM schemes 



We impose the momentum subtraction scheme 



Ai(/i) = l, 



(3.1) 



where 'Ai(/i)' stands for Ai evaluated at a specific kinematic point (e.g., symmetric or 
zero-momentum), with the momentum scale /i. The precise meaning of this will be clear 
when we discuss the MOM and MOM schemes. It then follows from ( 2.17| ) that 



g R (fi) = Z g l (p;a)g (a) = Z 2 (p; a)zl /2 (p,; a)g (a) Al ^ ^ ^ 



Z 2 (i,;a)Z 1 3 /2 (^a)go(a)Xr e (^a) 



(3.2) 



bare / 



go\i e is the quantity we calculate on the lattice. 

We will define two different renormalisation schemes, MOM and MOM. The 'asym- 
metric' MOM scheme is defined by setting the gluon momentum q 2 to zero. This differs 
from the MOM scheme defined in 0, and also from the MOM g scheme defined and com- 
puted to three-loop order in ||], where in both cases one of the quark momenta has been set 
to zero. The 'symmetric' MOM scheme is defined by the kinematics p = —k = —q/2 = s, 
so p 2 = k 2 = q 2 /4 = s 2 . The fully symmetric scheme where p 2 = q 2 = k 2 is impossible 
to implement on a finite lattice where the boundary conditions are different for fermions 
and gauge fields (antiperiodic and periodic in time respectively), which is why we are not 
considering it here. 
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In any scheme, the first step towards extracting the running coupling, which is propor- 
tional to Ai, is to eliminate those form factors with a different Dirac structure by tracing 
the vertex with 7^. With this in mind, we define the functions H^(p,q) as 



H^{p,q) = --lmTr 7 A1 A At (p, q) 



5o{Ai - (2p + g) 2 A 2 + [{p ■ q)q^ - q 2 p^](2p + q)^r 2 

- (q 2 - <il)n - [q ■ (2p + g) - q^p + ^]t 6 } , 



(3.3) 



where no sum over is implied. 

In the MOM scheme, all terms in (3.3) proportional to q or q^ disappear, and we are 
left with 

H„(p,q = 0) = g [x 1 (p 2 ,0,p 2 ) -4pl\ 2 (p 2 ,0,p 2 )) . (3.4) 

We can then eliminate A2 by imposing an appropriate kinematics: p^ = 0,p u 7^ for fj, ^ v. 
This defines Xf^fjj) = A^ 2 ,0, /1 2 ). 

The MOM renormalised coupling is then defined as 

^mom^) = Z2 ^) Z 3 1/2 (^oAi( / u 2 ,0,/, 2 ). (3.5) 



In the MOM scheme, the transverse-projected vertex gives us 

H T (p-2p) = g (4p 2 - 4p 2 j(-\A 1 (p 2 ,4p 2 ,p 2 ) - r 3 (p 2 , Ap 2 ,p 2 ] 



4p 2 

g {l- P ^)x[(p 2 ,4p 2 ,P 2 ) 



(3.6) 



A^ can then easily be extracted by 

X[( P 2 AP 2 ,P 2 ) = i^il>,-2p) ee h^p 2 ). (3.7) 



3 



Thus, we define the MOM running coupling as 

9mom(m) = Z 2 (^)zl /2 (n)g X' 1 (n 2 ,4fi 2 ,fj 2 ) . (3i 



4. Lattice formalism and computational details 

In this study, we use the Sheikholeslami-Wohlert fermion action, 

Ssw = Sw — i-rgoCsw ^2^2ij(x)a^F^(x)^(x) , (4.1) 

which is on-shell C(a)-improved (Sw is the Wilson action), along with an off-shell im- 
proved 2 quark propagator Si, given by fl8| 

Si(x,y) = (1 + b q ma)S (x, y) - 2ac' q 5(x - y) , (4.2) 

2 For full off-shell improvement, there should also be a gauge dependent improvement term. The absence 
of this term will give rise to errors, potentially of 0(g 2 a) |}4j]. We assume that this is a small effect compared 
with other systematic errors. Setting this term to zero is also consistent with mean- field improvement, which 
is what we are using in this paper. 
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where the 'unimproved' quark propagator So is simply derived from the inverse of the 
fermion matrix M: So(x,y) = (M (x, y; U)}. 

We define the lattice gluon field A^, which in the continuum limit becomes aA c ° nt , as 

% W 2r. , .1 (4-3) 



e 



UM - Ul{-qj) - ~ Tr (u,(q) - 



2i^o 

In order to reduce lattice artefacts, we employ a tree-level correction scheme, as dis- 
cussed at length in |l8|, The relevant tree-level expressions, and the definition of the 
tree-level corrected vertex form factors, are given in appendix |B|. 

All the results in this paper have been obtained with the Wilson gauge action at 



(3 = 6.0 on a 16 x 48 lattice. Using the hadronic radius tq [35, 36] to set the scale, this 
corresponds to a lattice spacing a -1 = 2.12 GeV. The gauge fields were generated with a 
Hybrid Over-Relaxed algorithm, with configurations separated by 800 sweeps. The quark 
propagators have been generated using a mean-field improved SW fermion action, with 
c sw = 1.479, for one value of k = 0.1370, or ma = 0.0579. Details of the computation 
are given in [37|. For the improvement coefficients b q and c' q of (fL^) we have used the 



mean-field values b q = 1.14, 2c' q = 0.57. 

The gauge fields have been fixed to Landau gauge, using a Fourier accelerated algorithm 



1 38] to deal with low- momentum modes. The Landau gauge condition has been achieved 



to an accuracy of yjr- u \ d^A^\ 2 < 10 12 . Further details of the gauge fixing are given 



m 



For the quark fields, we have used periodic boundary conditions in the spatial direc- 
tions and antiperiodic boundary conditions in the time direction. Hence, the available 
momentum values for an Nf x Nt lattice (with iVj, Nt even numbers and i = x,y, z) are 

K = ^(^"y) : ^ =1,2,- ■■ ,N; (4.4) 

2tt / 1 NA , \ 

Vt = ^[nt---^Y n t =1,2, • •• ,N t . (4.5) 



N t a V 2 2 

For the gluon fields, we have used periodic boundary conditions in all directions, and thus 
we have integer momentum values also in the time direction. The gluon tensor structure 



was studied in [13|. Our data correspond to the 'small lattice' in that paper. There it was 
found that 

T^(0)~diag(l,l,l,l/3) r(0)«y, (4.6) 



where T^ u and T are defined in ( |2.9|) . Significant deviations from the infinite- volume 
form were also found for the lowest one or two momentum points used for the symmetric 
kinematics. We have explicitly adjusted these points to account for this. For all other 
momentum combinations we will be studying here, the deviation of T^ y {q) from Pj^(q) = 
(Q(q)) were found to be negligible. 

5. Determination of Z 2 and Z 3 
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Figure 3: as a function of the renormalisation scale /ia, for Sq (left) and for Si (right), using 



the tree- level correction defined in (B.1C), and without any momentum cuts . 



In order to determine the quark field 
renormalisation constant Z2, we have used 
the tree- level corrected function Z(p), de- 
fined in ( |B.10| ). The results, for both So 
(for which Z^\p) = 1) and Sj are shown 
in figure [3|. 

In order to determine the gluon field 
renormalisation constant Z3, we use the 
simple tree-level correction procedure that 
was applied in [jDJ] . The gluon propagator 
D{q 2 ) is expressed in terms of the 'lattice 
momentum' Q [see QB.2| )], and a 'cylinder 
cut' is applied to select momenta near the 
4-diagonal. This is shown as a function of 
Qa in figure ||. We have fitted the gluon 
propagator to the phenomenological curve (Model A) of 1 13] 




Figure 4: Z3 as a function of the renormalisation 
scale ^a, for 125 configurations. The line denotes 



the best fit to the functional form (5.1) 



D(Q 2 ) 

where dp 
are 



AM 



2n 



+ 



{Q 2 + M 2 ) l + a Q 2 + M 2 



\n({Q 2 + M 2 ){Q- 2 + M- 2 )) 



-d D 



, (5-1) 



13/22 is the gluon anomalous dimension. The best estimates for the parameters 



2.02; 



A = 10.7; 



M = 0.534; 



2.17. 



(5.2) 



It should be emphasised that this fit is only performed to facilitate the computation of the 
running coupling, and no physical significance should be attached to the phenomenological 
parameters quoted. 
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6. Ai and the running coupling 



6.1 Asymmetric scheme 

We have calculated the proper vertex in the asymmetric scheme using both the unimproved 
quark propagator So and the improved propagator Si. 

We have evaluated Ai by first calculating H{{p,q = 0)(i = 1,2,3) for different values 
of p (with pi = 0), and then used invariance under the (hyper-)cubic group to perform a 
Z% average over i and equivalent values of p^ (as well as positive and negative p^ values). 
But first, we want to verify that the cubic invariance really holds. As figure || shows, all 
the three spatial components of the (uncorrected) vertex do indeed behave in the same 
fashion, within errors. The discrepancies of the order 2a can be put down to correlations 
between data at different momenta, combined with insufficient statistics. 

When p^ 7^ 0, H^p^q) also receives a contribution from A2. This means that we 
should not expect H4 to behave similarly to the other three components, since p^ = pt 
is necessarily non-zero. The lower panel of figure || confirms this — although part of the 
difference may also be due to finite volume effects affecting spatial and time directions 
differently. The form factor A2 will be studied in a forthcoming paper. 





Figure 5: Top: ifi (p, g) , ffj (p, 3) and H 3 (p 7 q) for q = 0, p = (0,p t ), as a function of p t a; using 
So (left) and Si (right). Bottom:_ff4(p, q) for q = 0, p — (0,pt), as a function of pta, for 83 
configurations using So (left) and for 100 configurations using Si (right). Note the different vertical 
scales for the upper and lower panels. 
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1 - 



0.5 1 1.5 2 2.5 3 

pa 

Figure 6: The unrenormaliscd form factor Xi(p 2 ,0,p 2 ) as a function of \pa\, with equivalent 
momenta averaged. The form factor taken from the improved propagator Sj is shown both before 
and after tree-level correction. After tree-level correction, the lattice data for Si lie on a single 
smooth curve. 



In figure || we show the unrenormalised form factor Ai, obtained by averaging all 
Hi(q = 0,pi = 0) over equivalent momenta and directions, as a function of \pa\. As the 
figure shows, this is a well-defined function of p (within the statistical errors) for pa < 1, 
both when So and Si is used. For larger values of pa, however, Ai extracted using Si 
developes significant ambiguities and a big 'bump' around pa = 1.7. This is due to the 
tree- level behaviour given in ( |B.15| ). Comparing with figure 1 of |^], we see that it is 
indeed approximately the inverse of the tree-level quark propagator. As expected, the 
irregular behaviour disappears after tree-level correction, and all the data lie on a single 
smooth curve. For pa > 1 this curve coincides with the data for So. 

We clearly see a substantial infrared enhancement of Ai, in accordance with the expec- 
tations from studies of the gluon and ghost propagators. Part of this must be an 'abelian' 
enhancement given by the QED expression (2.22) together with the infrared suppression 
of the quark propagator. We can determine the additional, 'non-abelian' enhancement 
by plotting the product of Ai and the quark propagator form factor Z(p), which in QED 
would be a momentum-independent constant. This is shown in figure [/], using the improved 
propagator Si. As we can see, a significant enhancement over and above the 'abelian' one 
remains, although our lattice is too small to allow us to draw any further quantitative 
conclusions. 

Using the values for Z2 and Z3 in section [|, we obtain <7jJqjj(^)) which is shown in 
figure ||, or, equivalently, a jfoM^' snown m fig ure H The coupling appears to reach 
a peak at about 1 GeV, below which it drops towards zero. However, caution is clearly 
warranted: the two or three lowest momentum points where this effect can be observed 
may well contain substantial finite volume effects (indeed, this was the case with the gluon 
propagator on the same lattice [fl3||), which only a simulation on a larger volume can resolve. 
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Figure 7: The unrenormalised form factor Xi(p 2 ,Q,p 2 ) multiplied by the quark renormalisation 
function Z(p), using the improved propagator Si, as a function of p. 
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Thus, at this point, we are not able to tell whether this may be a finite volume artefact, 
an artefact of the MOM scheme, or a real physical effect. Similarly, the fact that the peak 
value a ^j^ ~ 0.8 is very close to typical values for the frozen coupling extracted from 
phenomenology ||, may be suggestive, but nothing more. 

Turning now to the ultraviolet behaviour, we attempt to parametrise the leading non- 
perturbative and quark mass effects by fitting the results to the formula 

«(^-^=(l + ^)a 21oop (^), (6.1) 
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for [i > Pmin, where a 21oop is the two- loop running coupling, 

4vra 21oop (^) = \ , (6.2) 

y> 6 ln(/x 2 /A 2 ) + ^lnln(^ 2 /A 2 ) V ' 

with 60 = 11/16-7T 2 and b\ = 102/(167r 2 ) 2 the leading coefficients of the /3-function, for 
varying values of p m in- The results are shown in table ||. In all cases, the numbers obtained 
using So are almost identical to those obtained using Si, so we only report the latter. We 
may also, if we ignore the power corrections (i.e., set c = 0), compute A^^j directly 
according to the inverse of (16. 



A = ^e (b g 2 (fi)) 2b o . (6.3) 



The results of this are shown in figure 10. The numbers obtained by fitting this to a 
constant above p m %n are also reported in table |l|. These numbers are consistent with the 
result of fitting a(/x) to ( |6.2] ). We may also repeat this procedure after absorbing the power 
correction ( |6.1| ) in our definition of gn, using the value for c from our fit. The result of this 
is also shown in figure 1C and reported in table |l[ 



In lattice studies of momentum-space quantities, the momentum variable is to some 
extent arbitrary. We may choose any of the variables p, K(p), Q(p), K(p) or any other 
variable as long as it approaches p in the infrared and in the continuum limit, i.e. for pa C 1. 
If the continuum tree- level form of the quantity is momentum-dependent, we may use this 
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Pmin (GeV) 



11 



A (MeV) c(GeV 2 ) A (MeV) A r 



2.00 


155 


coo +63 


i.4±?:f 


oon+192 
•JO^-232 


425ifJ 


2.25 


149 


567ii 6 3 9 5 


i-i±f;l 


407^ 


454±gg 


2.50 


139 


555^141 


o.eifl 


445l|?| 


497±i| 


2.75 


125 


544+^ 


o.3±l;l 


/i«i+ 275 

461_304 


515±|i? 


3.00 


112 


534±}g 


-o.5±i;l 


O^-J-372 




3.25 


101 


5283?$ 


1 7+ 4 - 2 
-1. (-2A 


C1 o+366 
613_4i 


686_ 456 


3.50 


85 


536^182 


O +4.2 
-Z.Z_2.4 


662111 


7zLn +434 



Table 1: Fit parameters, using different momentum ranges. p m i n denotes the lower end of the fit 



window; the maximum in all cases being the maximum total momentum 5.75 GeV. n is the number 
of momentum points used in the fit. A is the value obtained for A^-^ without power correction, 
using (|6.3|), while A r is the value obtained by using the fitted value for c to extract a 21oop and 



feeding this into fl6.3|) 
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u, (GeV) 

(/Lt) (GeV) as a function of /i (GeV). Left: without power correction, j , . 
■ correction of (6.1) from a fit to /i > 2.0 GeV. The lines indicate the preferred 
value for A^^, with a 67% confidence interval. 



Figure 10: A 



= "MOM" ' \ - " ■ ' 

including the power correction 

^MOM 



Right: 



to guide our choice of variable [14]; however, when it is not, the choice remains largely 
arbitrary |pl] ]. Since the tree-level continuum vertex is momentum- indepentent, this is the 
situation we find ourselves in here. In order to quantify the resulting ambiguity, we have, in 
addition to the 'naive' momentum p, performed fits using K(p) = \/^2 u sin 2 (ap^) /a, which 



appears in the tree-level lattice vertex ( |B.12 ), as well as K z (p) = K(p)/Z^°\p), which is the 
momentum variable that makes the tree-level quark propagator take its continuum form. 
The use of this variable may be justified because the correction factor (p) appears also 



in the tree-level vertex, and also from the Ball-Chiu relation ( 2.22 ). The results of the fits 
are given in tables [2] and ||. 
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K m in (GeV) 


n 


A (MeV) 


c(GeV 2 ) 


A (MeV) 


A r (MeV) 


2.00 
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436+ 1 5 o 9 g 


o.4i?| 


343ti? 


374l|l 


2.25 
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423ii 


-o.2±i:§ 


404lli 
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423_ + 1 9 3 °o 
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78 


421 + 1 9 4 8 1 


-o.8±?:I 
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511^8 


3.00 


46 


425i|i 


-3.o±?;g 
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805lil 


3.25 


24 


455±iS 


-z.8_2.4 
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780lH 


3.50 


9 


470ii! 


O.Otfi 


42ll|li 
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Table 2: As table using K(jp) as our momentum variable. The maximum available momentum 
here is 3.70 GeV. 



K z>min (GeV) 


n 


A (MeV) 


c(GeV 2 ) 


A (MeV) 


A r (MeV) 


2.00 


155 


C1 r +65 

615_i 34 


1 1+ 2 - 


442li 2 


483l|l 


2.25 


150 


605 + 1 7 4 6 5 


i-i±2:i 


443112? 


484l|| 


2.50 


145 


596+i 8 4 3 9 


0.8! 2 ] 


461111 


503l|l! 


2.75 
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591_ + 1 9 6 4 4 


-1 0+2.6 


429li? 


.^+254 
466_283 


3.00 
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^0+106 
0OO_i72 


1 q+ 3 - 2 
-L-y-2.5 


39511? 


430111 


3.25 


112 


5781$ 


o o+4-3 
^••^-2.8 


37711? 


41 n +268 


3.50 


102 


571±ii2 


2.4±|5 


369li? 


4021271 



Table 3: As table [l], using K z (p) as our momentum variable. The maximum available momentum 
here is 5.05 GeV. 

From figure ^ and the right-hand panel of figure [l0| it would appear that the data are 
very well represented by a power-corrected two- loop running coupling as in (|0|), all the 
way down to 1.5 GeV if not lower. However, a glance at tables [l]-^ reveals several problems 
with this. 

Firstly, the fits are nowhere near stable. As the starting point for the fits goes from 2.5 
to 3.5 GeV, the best value for A increases by 50% when using p as our momentum variable, 
and the power correction goes from positive to negative. Secondly, the 'refitted' value for 
A, although always perfectly consistent with that obtained from Q6.1D , is consistently about 
10% higher. 

Thirdly, the fit values depend critically on which momentum value is used. To some 
extent this is simply because the values of p, K(p), and K z (p) may be very different 
when pa 3> 1, so different data are included in the fits. This is reflected in the different 
number of points for the same numerical value of the starting momentum. However, it 
also reflects a deeper ambiguity due to the finite lattice spacing. I.e., although there is no 
violation of 0(4) symmetry or other obvious signs of lattice artefacts in our data, and the 
near-perfect agreement between the So and Si results may be taken as an indication that 
lattice spacing errors are very small, those errors that do persist make a determination of 
a sensitive quantity such as A^^j prone to large uncertainties. It should be noted that 
we observe considerable anisotropy in the high-momentum region when gn is plotted as a 
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Figure 11: The unrenormalised, tree-level corrected form factor A'^p 2 , 4p 2 ,p 2 ) as a function of the 
quark momentum \pa\, from 498 configurations. 



function of K(p) or K z (p) as opposed to p, indicating that these are not the appropriate 
momentum variables in this case. Only by repeating the simulation at a smaller lattice 
spacing can this issue be properly resolved, however. 

Taking all this into account, we take as our best estimate for A the average of all 
the fits starting from 3.0 GeV (both with and without the power correction). This gives 
A.MOM = 530i|2o ± 100 ± 50 MeV, where the first set of errors are statistical, the second are 
due to the ambiguities in the choice of momentum variable, and the third is the intrinsic 
10% systematic uncertainty in the lattice spacing in the quenched approximation. 

6.2 Symmetric scheme 

At the symmetric point 2p + q = we use only the improved propagator Si, and thus the 
form factor A^ receives substantial tree-level correction according to ( |B.25| )-( |B.28| ). The 



tree- level corrected result is shown in figure 11. In order to reduce the statistical noise, we 



have averaged data for nearby momenta, within Apa < 0.05. We see that the data ara still 
considerably more noisy than for the asymmetric Ai of figure ||, but exhibit qualitatively 
the same behaviour. It appears that A'iO 2 ,^ 2 ,p 2 ) is more strongly infrared enhanced than 
Ai(p 2 ,0,p 2 ), but the noise makes it difficult to draw any definite conclusion. 

The MOM running coupling g MOM (^) is shown as a function of the renormalisation 



scale /x in figure [L2| The most obvious difference from the MOM coupling of figure |§| is 



that the noise is far worse, and we are not able to get any signal for A from these data. 



7. Matching to MS 

The relation between the scale parameters in two renormalisation schemes A and B is given 

by § 

Cab 



A^ = 6XP 



2&o 



(7.1) 



-17- 



2.5 



1.5 



0.5 



9_ <P 
o 



(1 (GeV) 



Figure 12: g MOM (n) as a function of /i (GeV). 



where Cab is the one- loop coefficient in the expansion of the coupling g 2 B in terms of g\. 

The complete one-loop expressions for the relevant form factors in the two kinematics 
we are studying, are given in appendix [C]. Here we only reproduce the results for the 
running coupling in Landau gauge, for Nf = 0. The MOM coupling is given by 



,151 3 m 2 9 ,/ m 2 



4 /i 2 

2 
2 



+ — In 1 + - 



4 v /x- 

2 * -4 25 m 2 



3 + 12 M 2 



16vr 2 



+ 0(g 4 



• (7-2) 



However, at asymptotically large momenta, where one-loop perturbation theory becomes 
valid, the corrections due to the mass term can be ignored. From ( |7.1| ) we thus find that 

Ar 



MOM 

A 



exp 



Using our 'best value' for A 



MOM 
N f =0 



MS 

we obtain 



151 
264 



1.77. 



A. 



MS 



300±i|8 ± 55 ± 30MeV. 



(7.3) 



(7.4) 



These numbers are above those obtained by other methods |q, 29], which yield a 'world 



average' of A. 



Nf=0 
MS 



240(10) MeV. With our large statistical and systematic uncertainties, 



our value is however fully consistent with the 'world average' 



The MOM coupling in the massless limit is 

SmomM =2Ms(/ i ) j 1+ (j ln2 + 
which gives for Aj^j 



793\ gMg(>) 
72 J 16vr 2 



0(g 4 



A 



MOM 
^MS 



oxi)i 8 1 n2/9 + 793/36A 28Q 



(7.5) 



(7.6) 
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If, instead, we renormalise the vertex at the gluon momentum, we find 

% = exp( 891n2 ^ + 793 / 36 U3.72. (7.7) 



"■MS 



8. Discussion and outlook 

We have studied the quark-gluon vertex in the Landau gauge, in quenched QCD with 0(a)- 
improved Wilson fermions, at two different kinematical points: an 'asymmetric' point, 
where the gluon momentum q is zero, and a 'symmetric' point, where q = —2p, in other 
words the incoming quark has equal and opposite momentum to the outgoing quark. 

We have focused on the form factor Ai, which is proportional to the running coupling 
and, in the decomposition given by ( A.6| ), ( |A.8j ), is the only form factor that is expected to 



be ultraviolet divergent. At the symmetric point, we are unable to study this form factor 
directly, and examine instead the linear combination A^ = Ai — q 2 T%. We observe that 
in both kinematics, Ai(A^) is substantially enhanced in the infrared. At the asymmetric 
point, this enhancement is significantly stronger than that expected in QED due to the 
well-established enhancement of the quark propagator form factor A{p). At the asymmetric 
point, no such direct comparison with QED is possible due to the admixture of T3, which 
is left unconstrained by the Ward-Takahashi (or Slavnov-Taylor) identity. However, the 
qualitative picture is the same. 

The lattice volume in this study is relatively small (a spatial length of ~1.5 fm and 
a total volume of 15-16 fm 4 ), so the infrared behaviour may well be contaminated by 
substantial finite volume effects. Although we have explicitly accounted for the large finite- 
volume effects appearing in the tensor structure of the gluon propagator, we have no 
guarantee that there are not substantial residual effects that, at the asymmetric point, 
could play an important role for all momenta. However, excellent agreement with results for 
Af-jg obtained by other methods have been obtained from the three-gluon vertex in a MOM 
scheme on symmetric lattices [j2S[j, and there appears to be no reason why the situation 
should be much worse in our case. The qualitative similarity between the symmetric and 
asymmetric point might also be taken as an indication that finite volume effects, although 
possibly sizeable, do not dominate. In any case, it would be desirable to perform the 
simulation on a larger lattice in order to have a better resolution of the momentum in 
directions other than the time direction. This is also the only way we would be able to 
settle the issue of whether the running coupling is frozen, or possibly goes to zero. 

We have used the results for Ai at the asymmetric point to determine the running 
coupling a s in a zero-momentum (MOM) renormalisation scheme, and obtained from 
this a nonperturbative estimate of A^L °. Our main results are for the strong cou- 
pling a^C?(2GeV) = 0.36(4); c^ =0 (2GeV) = 0.28(3), and for the QCD scale A^ =0 = 

300^180 ± 55 ± 30 MeV, where the first set of errors are statistical, the second due to ambi- 
guities in defining the momentum variable, and the third due to uncertainty in the lattice 
spacing. This is consistent with, although slightly higher than other estimates for A^g. 
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Although the excellent agreement between the results using the 'unimproved' propa- 
gator So and the 'improved' propagator Si indicate that our tree-level correction scheme 
has successfully accounted for the large high- momentum lattice artefacts, and that resid- 
ual 0{a) errors are not a significant factor, the need for large tree-level corrections still 
implies some uncertainty about the results, at least in the intermediate momentum regime. 
A fermion action which is more well-behaved at high momenta, such as overlap fermions, 
would be a great improvement. 

The main source of systematic uncertainty, and of possible discrepancies between our 
result for A— £ ° and those of other determinations, is that we have not been able to access 

MS ' 

sufficiently high momenta, where two-loop scaling should be valid, nor have we taken into 
account higher-order perturbative effects, which should extend the range of validity for 
the perturbative matching. Experience from the 3-gluon vertex [29| suggests that both 
a large momentum window and 3-, perhaps 4-loop running of the /3-function are needed 
to obtain reliable results. This requires simulations at smaller lattice spacings, as well as 
a two-loop calculation of the Ai form factor in the relevant kinematical limit. Both are 
computationally very expensive. 

As we mentioned in the introduction, two-loop calculations have already been per- 
formed in both an asymmetric 0] and a symmetric || kinematics. Neither is, however, the 
kinematics we are employing here. 

In the MOM (symmetric) kinematics, we have been unable to get a reasonable signal 
for the running coupling. The main reason for this is statistical noise, but the need for 
large tree-level corrections is clearly also a significant factor. For this reason it would 
be essential, if we were to attempt a more accurate determination of the vertex in this 
kinematics, to choose a fermion discretisation which is not afflicted by such problems. 

In a precision study, the quark mass must also be handled carefully. Here, we have 
merely included the quark mass in the overall power correction, which has been determined 
numerically. An obvious next step would be to study the vertex at a second quark mass. It 
would be an advantage, also for this purpose, to use a fermion action which respects chiral 
symmetry, such as overlap fermions, or a remnant thereof, such as staggered fermions. 

The large numerical uncertainties have prevented us from obtaining any reliable esti- 
mate of the power correction. An alternative approach would be to calculate analytically 
the size of the power corrections from the condensates involved, of which the dominant is 
expected to be the chiral condensate, using the available estimates for the values of the 
condensates. 

The running coupling extracted in the MOM scheme reaches a maximum at 0.8 GeV. 
A similar result was found for the three-gluon vertex in an analogous MOM scheme |§,f2j|. 
This would correspond to a zero in the /3-function at the maximum coupling, with double 
values below that. It has been suggested [40] that this can be related to infrared singu- 
larities in the ghost self-energy. Such singularities should not affect symmetric momentum 
subtraction schemes. Our results in the MOM scheme can neither confirm nor refute this 
conjecture. 

In order to resolve this issue, and to pin down the low- and intermediate-momentum 
behaviour of the quark-gluon vertex, simulations on larger, and possibly coarser lattices 
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are necessary. This is an orthogonal line of inquiry to that needed to determine the running 
coupling along with the power corrections, which requires much finer, but not larger lattices. 

Work is currently in progress to determine the form factors A2 and A3 at the asymmetric 
(q = 0) point. These form factors both vanish at tree level in the continuum, but must 
be non-zero nonperturbatively in order to fulfil the Slavnov-Taylor identity. A complete 
determination of all form factors would be a natural next step. This is, however, not 
possible in Landau gauge because of the transversality condition. For this reason, and also 
because the gauge dependence of the vertex is in itself of theoretical importance, it would 
be of great interest to study the vertex in a generic covariant gauge |4l], ^2]. This would 
also allow calculations in the unmodified symmetric MOM scheme, which might have some 
advantages over the modified scheme we have used here. 

At present, there is no known method to reliably assess the effect of Gribov copies. 
All numerical methods founder on the fact that as the physical volume increases, the 
number of Gribov copies also increases, and it becomes impossible to ascertain that one 
has found either the absolute maximum of the gauge fixing functional, or any other unique 
representative. Choosing a gauge without Gribov copies, such as the Laplacian gauge [43] 
or axial gauges, does not solve the problem, since results in one gauge tell us nothing about 
the effect of Gribov copies in a different gauge. 

On a practical level, attempting to select, however imperfectly, the absolute maximum, 
using e.g. 'brute force' |24|], simulated annealing |lq| , or smeared gauge fixing [44] is in 
principle worthwhile. At present, however, we would expect any signal showing a difference 
between the naive (minimal) Landau gauge and the fundamental modular domain to be 
swamped by statistical noise for three-point functions such as the quark-gluon vertex. 
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A. Tensor decomposition of the vertex 

The Lorentz structure of the vertex in the continuum consists of 12 independent vectors 
and can be written as 

12 

= -ig T„ = -ig Yj > ( AA ) 

i=l 
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where 
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(A.2) 



It is useful to divide the vertex into a 'Slavnov-Taylor' (non-transverse) part and a trans- 
verse part, as is commonly done in QED: 

A>, q) = A(f T )(p, q) + A^\p, q) . (A.3) 



The ST part is that part of the vertex that saturates the Slavnov-Taylor identity ( |1.2j ) 
and contains no kinematical singularities. It is often, misleadingly, referred to as the 
'longitudinal' part, although it also contains a transverse component. 



We will make use of the QED decomposition ||, £6], 27] of the fermion-gauge-boson 
vertex function, which is usually given in Minkowski space. We wish to write the euclidean- 
space equivalent, in such a way that all the scalar form factors Aj and Tj are the same as 
in Minkowski space. The usual procedure is to apply the Wick rotation (jjq — ► ip4,Pi — ► 
—Piilo ~^ 74; li ~^ but since the vertex is a four-vector, this is not a linear trans- 
formation in our case. Our prescription is to create a Lorentz scalar by contracting the 
vertex with 7^ and require that the Wick-rotated Minkowski result is identical to what 
we obtain by performing this operation in euclidean space. In particular, any euclidean 
scalar function should be equal to the Minkowskian scalar function sampled at spacelike 
(i.e. negative) momenta: 

f E ( P 2 ,q 2 , k 2 ) = f M (-p 2 , -q 2 , -k 2 ) , (A.4) 

where / can be t, or Aj. 

Following this prescription, we can write the ST part as [with k = p + q] 

4 

Af%, q) = ~ig Hp\ Q 2 , k 2 )L itfl (p, q) , (A.5) 
i=l 

where the euclidean-space functions Fj jAt are given by 

Li )fl = 7m ; L 2 ^ = -{20+ i) (2p + q)^ ; 

Lz,n = -i(2p + q)n; L 4jjU = -ia^ (2p + q) v . (A. 6) 

For the purely transverse part PS T \ we will use the decomposition of [^7|, which 
differs slightly from the one of H [2f|. This decomposition is preferable because it is free 
of kinematical singularities in all covariant gauges. Moreover, as we shall see, the relation 
between the ST and purely transverse parts of the vertex becomes more transparent in this 

(T) 

basis. The purely transverse part of the vertex is specified by q^A^ (p, q) = and satisfies 
A/t '{p, 0) = 0, and we write 

8 

A^fa q) = -ig^^P 2 ' I*' k2 ) T ^ q) . ( A -7) 
i=i 
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where the euclidean-space functions Tj are given by 

Ti tfl = i [p^q 2 - q^p ■ q)] ; 

T20 = [p^q 2 ~ q„ip ■ q)] (2 jf+ $) ; 

T 3,^ - ? 2 7m ; 

= -i [g 2 cr MZ ,(2p + g) v + 2q (1 a v \p v q > \ ; 
r 5) „ = -i<V?„ ; (A - 8) 
T 6 ^ = q-(2p + q)^- 4{2p + q)^; 

T 7,n = -^q ■ (2p + q) [(2 jf+ - (2p + - i(2p + q) lM a uX p v q\ ; 

^ = -l^xPvQx- fan+ iPn ■ 

Charge conjugation symmetry dictates that all the Aj's and Tj's are even with respect to 
interchanges of p 2 and k 2 (or (p + q) 2 ), except for A4,T4 and tq, which are odd. 

In this decomposition, the transverse projection of the ST part of the vertex is given 

by 



P^u{q)Li, u (p,q) = 


-■^T 3tli (p, q); 


(A.9) 


Pfj,u(q)L 2 ,u(p,q) = 


2 

—^T 2 ^(p, q); 


(A.10) 


P^,u(q)L3,u(p,q) = 


2 , ^ 

—^ T iAP> <?); 


(A.ll) 


Pfj,u(q)L4,u(p,q) = 


j2-T 4 ^(p, q). 


(A.12) 



Thus, we will define the following modified form factors, which appear in the transverse- 
projected vertex: 

A; = Ai - q 2 r 3 ; A' 2 = A 2 - ^r 2 ; (A.13) 
q 2 

•Ms = ^3 _ y r i ; A4 = A 4 + ^ r 4 . 
B. Tree-level lattice expressions 

We define and use the following momentum variables, which may be used to bring the 
lattice tree-level expressions into a more continuum-like form, 

K nip) = ^sin^a), (B.l) 



2 / 

Qui?) = ~ sin(p^a/2) = — yjl - cos(p M a) , (B.2) 

K^p) = ^(2p) = i- sin(2^a) , (B.3) 

CjiGO = cos(p M a) . (B.4) 
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At tree level, the dimensionless momentum-space propagator So(p) is identical to the 
free Wilson propagator, 



-ia JC(p) + ma+ \a 2 Q 2 (p) 



a 2 K 2 (p) + (ma + \a 2 Q 2 (p)) 
The tree- level form of the C(a)-improved propagator Si is given by 



Sf> (p) = (1 + b q am)S^ (p) - 2ac' q = — 

ia fC(p) + amZm (p) 



where 



1 + b q am 
DM 



K 2 (p) + {m + -Q 2 {p)Y 



with 



(B.5) 



(B.6) 



(B.7) 



D/(p) = (1 + b q am) 2 K 2 (p) + B 2 (p) , 

Bj(p) = (1 + 6 9 am)(m + 1q 2 ( P )) - 2a< [# 2 (p) + (m + \q 2 ( P ))' 



(B.8) 
(B.9) 



The tree-level functions Z^ ) and Zm' (which is given in [JH]] and will not be reproduced 
here) give rise to very large finite-a effects at intermediate and large momenta. To reduce 
these lattice artefacts, and bring the high- momentum behaviour of the quark propagator 
into contact with the continuum perturbative behaviour, we employ the tree-level correction 



scheme defined in [18, [LS|], where the quark propagator is written as 

1 



s-\ P ) 



Z{p)Z(°){p) 



ia JC(p) + aM h (p)Z { + ) (p) + oAM ( -'(p) 



(B.10) 



with the functions Z^ , Zffi and AM ( ~' denoting the tree-level behaviour. We call the 
functions Z(p) and M h (p) the tree-level corrected quark form factors. Here, we are only 
interested in Zip), which can be related to the quark field renormalisation constant Z%, so 
we can ignore the mass correction functions Zm\ AM (_) which are defined in | |19| |. 
At tree level, the Landau gauge gluon propagator becomes 



D^q) = P$(q)D(Q 2 ) = (s^ 



Q»(q)Qv(q) 



l 



Q 2 (q) J Q 2 (q) ' 

In this notation, the gluon propagator requires no further tree-level correction. 



(B.ll) 



The tree-level lattice vertex using the 'unimproved' propagator So is [45, 46 



<» M = 7, cos a -^P^ - isiu a ^±- 



aqti 



fiu cos sin aq v 



(B.12) 
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The constant term c' in Sj does not contribute to the unamputated vertex in (]2.4|), 
since {A^) = 0. Thus, the improved vertex at tree level is given by 

A?JW = (l + b q am)sP(p)- l sl?\p)A a ^\p, q )sl?\p + q )SP(p + q )- 1 . (B.13) 

The full expression is very complicated, but it simplifies greatly for the two cases (symmetric 
and asymmetric) in which we are interested. 

B.l Asymmetric kinematics 

In this case the gluon momentum q = 0, while the quark momentum is 'orthogonal' to the 
vertex, i.e. the //-component p^ of the quark momentum is zero. Then the tree-level lattice 
vertex ( |B.12|) reduces to 

At%,0) = -; 5 ot a 7M- (B.14) 
Making use of this, along with the unimproved flB.5j) and improved (|B.6|) propagators, the 



improved vertex (B.13) becomes 



O,°)W=0 = (l + b q am)sf\ P )- 1 S^(p)AZ ) (p,0)sS ,) (p)sP(p)- 1 

(B.15) 

= -ig t a lll /Z^(p). 
Thus, at this point, the tree-level corrected vertex may be defined according to 

A?>'°)k=o = -^ a ^^y5oAi(p 2 ,0,p 2 ) 7At . (B.16) 

As previously mentioned, in the Landau gauge we calculate the transverse- projected 
vertex ( |2.12 ). In the asymmetric case, this becomes 



P^(q)K(p,q = 0) = VAf'feO) + AjPfoO) = 4 ST) (p,0) , (B.17) 



since (p,0) = and P^(0) = S^ u as discussed on p. ||. The lattice, finite- volume 
version of this is 

T^(q)A v (p,q = 0) | PM=0 = T^(0)A(f T )(p,0) |^ =Q = -i^^-g^^, 0,p 2 ) 7/1 . (B.18) 
B.2 Symmetric kinematics 

The second case we consider is the (symmetric) case where 2p + q = 0, i.e. p = —koi 
equivalently q = —2p. In this limit, the tree- level vertex ( p.!2| ) becomes (for ease of 
notation we will here set the lattice spacing a = 1) 

A^J? (p, -2p) = -igvVQ (p, -2p) = -ig t a ^ + ic sw £ a^C^(p)k v (p) \ , (B.19) 

and repeating the same procedure as in the asymmetric case, the improved vertex ( |B.13| ) 
takes the form 

A^(p, -2p) = (1 + b q m)SP(p)- 1 sl°\p)A a ^(p, -2p)sS 0) (-p)5f ) (-p)- 1 . (B.20) 
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This gives us 
D 2 



i "! T?l(p, -2p) = B v - A 2 V K 2 + 2A v B v c sw (K-K)C^ 7/i + 2A 2 V JiK^ 

- 2A yJ Bv/c sw fCKp - 2iA v B v ^ a^i^ + ic sw {A 2 v K 2 + ^)C„ ^ a^k v 

V V 

- 2ic sw A 2 v (K-K)K^ + 2ic sw A 2 v (K-K)C^ - 2ic sw A 2 v K^ £ a vX K v K x , (B.21) 

where we have written if = K(p),K = K(p), C = C(p) and 

A v = A v (p) = (i + 6 gm )( m + Ig 2 ( p )) -Bj(p), (B.22) 

B v = B v (p) = (1 + b q m)K 2 (p) + (m + ^Q 2 (p))B/Cp) . (B.23) 

If we concentrate on the part of this that becomes proportional to Ai and T3 in the contin- 
uum, the lattice expression may be decomposed as 

Vfl(p, -2p) = AfV - ±4°\K 2 ^- fCK,,) - a4 0) (K-KC^- fCK„) + ... (B.24) 

We can read off the tree- level form factors from ( |B.21| ), 

\f\p 2 ,Ap 2 ,P 2 ) = l -j^(A 2 y{p)K 2 (p)+By{p))=l/Z^{p)- (B.25) 

rj°V,4pV)= ^^A 2 v (p); (B.26) 



4 0) (P 2 AP 2 ,P 2 ) = - 1 -^^c sw A v (p)By(p) . (B.27) 



As in the asymmetric case, we compute the transverse-projected vertex (|2,12 ) in Landau 
gauge. With the decomposition ( |B.24| ), it reads (for sufficiently large q) 

rj(p,-2p) ^ I**(q)r v <p,q= -2p) = - Q ^f P) ) r„(p, -2p) 



K ^{p)K v { p) 
K 2 (p) 



M/K* - 4r 3 {K z ltl - JCKJ - 4f 3 (K-KC^- fCK^) + ... (B.28) 



This is the lattice equivalent of the projection (2.24). From this we obtain the transverse- 
projected, lattice equivalent of Q3.7|), 



hi(4p 2 ) = --^ImTV 7/1 Aj(p,-2p) 



3 50 (A!(p 2 , Ap 2 ,p 2 ) -AK 2 r 3 (p 2 , Ap 2 ,p 2 )) ( B - 2 9) 
- Ag [aK-K -K 2 - ^Q 2 K-k)Up 2 , Ap 2 ,p 2 ) = 3<7 A? at . 
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The tree-level corrected A' x is therefore 

A'^V, 4p 2 ,p 2 ) 



A'i(p 2 ,V,p 2 ) 



A'/°^(p 2 , 4p 2 ,p 2 ) 



MV)/<?o 



(B.30) 



>(o) 



4^ 2 ri^ 



4(4K-K -K 2 - \Q 2 K-K)ff ] 



3(A X - , 3 

where the p-dependence in the denominator on the last line is implicit. 



C. One-loop expressions 



In this section all the expressions will be given in Minkowski space. In the MS scheme, 
the one-loop contributions S^,IlW an d A^ to the quark and gluon self-energy and the 
vertex component Ai in the MOM kinematics are given by [0, 0] 



16^ 2 I S 



2 2 2 4 / 2 

m m — p m I p 
H 7T - m 7. 1 T In 1 

p z /i z p 4 \ m z 



m 2 / m 2 \, / p 2 
p \ p J \ mr 



(C.l) 



nW(p 2 ;M ) 



16vr 2 



"36-2^ + U-2J ln ^ 



21 



C A 



+ \T R N f 



1 / ^rn 2 \ , m 2 
3( 5 + 12 ^J" ln ^ 



+ 1 + 



2m z 



p- 



4m z 
p 2 



2\ 1/2 



hi 



4/7 r 



f 1 — 4m2 



1/2 



+ 1 



16vr 2 



m 



ec F + (3 + e) 



In 



(3+o + (i-0— 

2 2 

m — p 



+ 



ec F + (i-o 



4 / 2 

m / p 

— In 1 2 

p* \ ?7l z 



(C.2) 



(C.3) 



Setting p 2 = —/x 2 , this gives the following expression for the MOM running coupling 
9 mom^' in Landau S au S e (£ = 0) for N f = 0) 



1 + a« (V, 0, V; M ) - e? } (V; m) - (V; //) + 0( 5 4 ) 



1 + 



151 3 m 2 9 



- - In 1 + -TT 



24 4 fi 2 4 



m 



+ -=-ln(l + 



m- 



4 25 m 2 
3 + 12 M 2 " 



16vr 2 



• (C4) 
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In the MOM kinematics, the one-loop contributions to Ai and T3 are given by 
AiV,4 S V; M ) = 



16vr 2 



ho A +(CF+\c A )^+ 



m 



In 



2 2 

77T - S Z 



A* 



-^(i + O-^- 



, -4s 2 
In — — + 



I 1 



■(i + O- 



ln- 



A 4 



+ 



s" 1 + Am z s z + m 
s 2 (s 2 + m 2 ) 



m 



In(l--J) ; (C.5) 



1 



16vr 2 12s 2 



(2C F - CU) 



.(2C f --Ca)(2-0 



2 + 2£ + 



2 

m 2 (5 + 



s 2 — m 2 



l-^ln Vl 

s" 1 



4C F + C A (l-e) 
1 



m 
T2" 



/ 1 m 2 



+ ln^- 
A* 



+ < 4C F (1 + - ^(7 - 2£ - f) + (4C F - C A ) r °"~ 



m" 3m 

9~ ^ + 9 9 

s z s z — m z 



-C A 



m 



s 2 + m 2 



m 



In 



22 2 
m — s m 



A 4 



In 1 ^ 



s 2 + m 2 



(5 _4£-^) + —(! + £) 



2 (s 2 + m 2 ) 2 



[ s 2 ( 3 _ 6 £ _ £ 2 ) + m 2 (i 3 _ u£ - 3 £ 2 )] In • 



-4s 2 



In Landau gauge, for SU(3), we find in the massless limit that 

Af } (s 2 , 4s 2 , s 2 ; M ) = A< 1} (s 2 , 4s 2 , s 2 ; y) + 4s 2 r 3 (1) (s 2 , 4s 2 , s 2 ; y) 



<4,(a0 f2hl 4 9 -s 2 

+ - In 2 - - In — tt 

4 n 2 



16vr 2 V 36 9 



(C.7) 



Setting s 2 = — y 2 , analogously to the MOM scheme, we find the MOM running coupling 
at asymptotically large momenta to be 



5mom(^) — 0Ms(aO 



793 A. fMs(^) 



(C.8) 



A 4 



(C.6) 
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